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The dispersion relation of surface waves of a magnetic fluid 
in a magnetic field is studied experimentally. We verify the 
theoretically predicted existence of a non-monotonic disper- 
sion relation. In particular, we demonstrate the existence of 
two different wave numbers occuring at the same frequency 
in an annular geometry. 



I. INTRODUCTION 

The science of magnetic fluids is a multidisciplinary 
field, which encompasses physics, chemistry, mathemat- 
ics, engineering and medical sciences Q. One of its inter- 
esting features is the spontaneous formation of patterns 
under the influence of an external mag netic field g. If 
a static magnetic field is applied normal to the surface, 
the so called normal field- or Rosensweig-instability will 
occur at a critical value of the magnetisation leading to a 
deformation of the surface. Due to theoretical consider- 
ations |Q , the advent of this pattern forming instability 
is accompanied by a non-monotonic dispersion relation. 
This interesting behaviour is not a unique feature of mag- 
netic fluids, however. Similar behaviour is expected for 
polarizablc fluids under the influence of a normal elec- 
tric field ||. Both theories apply for an inviscous fluid 
in an infinite layer. Although experimental studies of ex- 
ternally induced surface waves have been published pl^] , 
an experimental verification of the non-monotonic disper- 
sion of surface waves is missing so far. We thus present 
an experimental attempt to verify the existence of a non- 
monotonic dispersion relation for surface waves in a mag- 
netic fluid. 



II. EXPERIMENTAL SETUP AND PROCEDURE 

The experimental setup is shown in Fig. |l|. We use 
the commercially available magnetic fluid EMG 909 (Fer- 
rofluidics) with the following properties: density p — 



1020 kgm 3 , surface tension a — 2.65 • 10 2 kgs 2 , ini- 
tial magnetic permeability (i = 1.8, magnetic saturation 
M s = 1.6 • 10 4 Am -1 , viscosity 77 = 6 • 10" 3 Nsm -2 . 
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FIG. 1. Experimental setup. 

The fluid is filled into a V-shaped || cylindrical teflon 
channel of 6 cm diameter as indicated in Fig. |2|. The 
channel has a depth of 6 mm and upper width of 5 mm. 
The motivation for using the V-shape is to minimize field- 
inhomogeneities at the boundaries of the channel. The 
upper part of the channel has a slope of 15°, which is 
close to the measured contact angle between EMG 909 
and teflon, in order to enforce a flat surface of the fluid. 

The channel is placed on the top of an electromagnet, 
with an iron core of 9 cm diameter and 15 cm height. 
The coil consists of 250 windings of copper wire with a 
diameter of 2 mm. A current of 4 A is then sufficient 
to produce the magnetic field used in these experiments. 
The two additional iron cylinders indicated in Fig. |l| on 
top of the teflon piece have a diameter 3.5 cm and height 
3 cm, and diameter 3 cm and height 3 cm, respectively. 
They serve to increase the homogeneity of the magnetic 
field near the surface of the fluid. The field is monitored 
by means of a hall probe (Koshava 3 Teslameter) located 
near the surface of the channel. 

In order to excite surface waves without mechanical 
contact with the magnetic fluid, we use an air jet created 
by a loudspeaker of diameter 10 cm, which is directed 
towards the surface of the fluid with the help of a noz- 
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zle as indicated in Fig. [|. The loudspeaker is driven by 
a sinusoidal current with computer-controlled amplitude 
and frequency. 
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FIG. 2. Cross-section of the V-shaped teflon channel. 

The shadowgraph method |7]] is used for detecting the 
surface waves. The light from a tungsten bulb, placed 50 
cm above the center of the channel, is reflected at the sur- 
face of the fluid and directed towards the CCD-camera 
with a glass plate acting as a semitransparent mirror as 
indicated in Fig. |l|. In order to avoid strong nonlinearities 
of the optical response, the lens of the camera focuses on 
a plane only a few millimeters above the surface of the 
fluid. A more precise determination of the focal plane 
is not necessary because no attempt is made to measure 
the amplitude of the surface deformation, as this is not 
considered as an important parameter for these experi- 
ments. 

The analysis of the images and the control of the ex- 
periment is done with a 66 MHz 80486-PC, equipped 
with a 512 x 512 8 bit frame grabber (Data Transla- 
tion DT2853), a programmable counter (8253) located 
on a multifunction I/O-board (Meilhaus ME-30), and a 
synthesizer-board (STAC ASP-ADW1). 

The CCD-camera (Phillips LDH 0600/00) works in the 
interlaced mode at 50 Hz using an exposure time of 40 
ms. The time between two horizontal synchronization 
pulses is Ths = 64 /is. The corresponding frequency 
of 1/Ths = 15625 Hz is divided by the programmable 
counter. The reduced frequency at the output of the 
counter acts as a pacer for the synthesizer-board, driving 
the loudspeaker via an amplifier (Euro Test LAB/S 115). 

The camera takes images at a fixed frequency of 25 Hz, 
while the driving frequency for the surface waves of inter- 
est reaches up to 30 Hz. The ensuing problem in resolving 
the dynamics is overcome by a software-realized strobo- 
scopic image acquisition, which is based on the computer- 
controlled synchronization between driving and sample 
frequency. 

The stroboscopic image acquisition works in the fol- 
lowing way: The period To of the driving oscillation can 
only be an integer multiple of the time Ths between hor- 



izontal synchronization pulses according to 

Td = NsbMcThs, 

where Mc is the multiplier value loaded in the counter, 
and Nsb is the number of data points used by the 
synthesizer-board to represent the sine wave. The im- 
ages can only be sampled with a time interval T$ = 
625 Ths — 40 ms, which is the time between two full 
frames. Consecutive images are taken at times tj = jTs, 
where j represents the number of the image. If one as- 
sumes the response of the surface to be periodic with the 
driving period Td, then the surface at time tj is in the 
same state as at time t'j < Td, according to: 

t'j{tj) = t^modTo = (jTs)modTb 

= ({625j)mod{N SB M c ))T H s- 

This equation describes the projection of the images 
taken at the times tj into the time interval [0,Td [• The 
ordering of these images according to increasing tj(tj) 
yields to the reconstruction of one period of the surface 
motion. The total measurement time Tm is the least 
common multiple of the two periods To and Tg. The 
number of different values of tj is Nj = Tm/Ts- The 
driving frequencies are chosen as a compromise between 
the desired frequency and sufficient temporary resolution 
corresponding to N] « 100. 

In order to measure the light intensity along the annu- 
lus by means of the frame camera, we divide the circle 
into 128 segments as schematically indicated in Fig. ||. 
The mean value of the intensities of those pixels which 
are located within one segment of the ring is calculated, 
thus yielding to 128 intensity measurements along the 
annulus. Wc consider only odd-numbered lines, which 
correspond to one half-frame of the camera working in 
the interlaced mode. The calculations can be performed 
with the video frequency of 25 frames per second by us- 
ing a look-up-table technique for the addressing of the 
pixels of a segment. 
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FIG. 3. The rectangular pixels of the frame camera are 
assigned to segments along an annulus. Only 16 segments 
are shown here for clarity; for the measurements we used 128 
segments. Note that the number of pixels per segment is not 
constant. 
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III. EXPERIMENTAL RESULTS 



FIG. 4. Space-time-plot for 5 = 0, and T D = 768 
Ths = 49.152 ms. 



and 10 below). In principal this experimental procedure 
would also allow to determine the attenuation of surface 
waves, provided that our optical setup responds in a lin- 
ear manner. 
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FIG. 5. Spectrum obtained from Fig. W. Only wave num- 
bers corresponding to a frequency of — 20.345 Hz are 
shown. 

Fig. ^ demonstrates the behavior of a surface wave 
which is driven at an excitation period To = 49.152 
ms in the absence of a magnetic field. The horizontal 
axis represents the 128 intensity values measured along 
the annulus. The driving occurs at 0° and 360°, respec- 
tively. In this case the stroboscopic algorithm yields 256 
measurements per period. These 256 lines are plotted 
on top of each other according to the sorting algorithm 
described above. The image displays two waves ema- 
nating from the source, travelling clockwise and coun- 
terclockwise along the channel. The small amplitude at 
the position of 180° results mainly from the attenuation 
and presumably partly from the interference of the two 
waves. In order to extract the wave numbers we per- 
form a two-dimensional FFT JsJ of this image. In this 
case this method is artifact-free because we have peri- 
odic boundary conditions both in space and in time. The 
FFT leads to a two dimensional complex field H(k,oj). 
The frequency of interest is ui\ = 2tt/To- In Fig. || we 
have plotted \H(k,oji)\ as a function of the 128 avail- 
able wave numbers. Negative wave numbers represent 
the wave travelling in clockwise direction, while positive 
wave numbers correspond to the ones travelling counter- 
clockwise. The spectrum clearly shows two peaks indicat- 
ing the two dominant waves. The finite size of the peaks 
is due to the attenuation of the waves. The wave number 
kp is interpolated between the peak value \H(k max i,uJi)\ 
and its largest neighbour \H (k maX 2, u>i)\ according to 

, k max i\H(k max i, u>i)\ + k maX 2\H(k maX 2,uJi)\ 

kp = . 

\H{k m axi,u\)\ + \H(k max2 ,uJi)\ 

Measuring kp for different values of the driving period 
Tp> then yields to the dispersion relation (see Figs. || 




FIG. 6. Fully developed peaks of a magnetic fluid obtained 
at B > B c in the annular channel. 

If a sufficiently high magnetic field perpendicular to the 
surface of the magnetic fluid is applied, the normal field 
instability occurs, leading to characteristic peaks shown 
in Fig. [| Because of problems with the aging of the 
magnetic fluid, the critical magnetic field B c for the onset 
of this instability drifts in a few weeks from B c — 0.009 
T to 0.011 T. We thus present all values of the magnetic 
field in units of the corresponding critical magnetic field. 

On the basis of theoretical considerations for the on- 
set of the Rosensweig-Instability Q one can show that 
non-monotonic dispersion relations are expected for val- 
ues of the magnetic field ranging between 0.93 B c and 
B c . This range is independent of the parameters of the 
fluid. An example for a non-trivial behavior of the sur- 
face waves in this regime is shown in Fig. Q for a value of 
the magnetic field of 0.93 B c . This space-time plot dis- 
plays a clearly more complicated behavior than Fig. ^. 
It seems to indicate that more than one wave number is 
excited by the driving period To = 196.608 ms. The cor- 
responding spectrum is shown in Fig. |[ The extraction 
of the dominant wave number is easily obtained, while 
peaks corresponding to higher wave numbers are much 
less pronounced and cannot be obtained from this spec- 
trum with reasonable accuracy. We have thus decided 
to extract only one wave number from the spectra. The 
results for counter-clockwise travelling waves are sum- 
marized in Fig. ||, and for the clockwise waves in Fig. [!(]. 
The ordinate of the plot corresponds to the driving fre- 
quency, while the abscissa indicates the wave number ob- 
tained from the spectra. The results for zero magnetic 
fields are shown as crosses. 



FIG. 7. Space-time-plot for B 
Ths = 196.608 ms. 



0.93 B c , and To = 3072 
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FIG. 8. Spectrum obtained from Fig. ^. Only wave num- 
bers corresponding to a frequency of T^ 1 — 5.09 Hz are 
shown. 

The solid lines in Figs. || and [l0| are fits of the function 
lo 2 = ak — bk 2 + cfc 3 , where a, b and c are fit parame- 
ters. The choice of this mathematical form for the fitting 
function is inspired by the dispersion relation obtained 
for an inviscid fluid and an infinite layer [Q. There the 
parameter a would have the value of g = 9.81 m/s 2 . The 
mean value of our fit parameter a is 9.96 m/s 2 . The same 

theory predicts b to be given by the term fiQji ^^ j;H 2 . 

The mean value of b/H 2 is 2.9 • KT 10 When using 

the permeability /i = 1.8 provided by the manufacturer 
of the ferrofluid, one expects 5.1 • 10~ 10 ji^r- The pa- 
rameter c is theoretically predicted tobeer//9 = 2.6-lCP 5 
m 3 /s 2 when using the values of the surface tension and 
the density given above. The measured mean value of 
c is 2.0 • 10~ 5 m 3 /s 2 . In spite of this fairly reasonable 
agreement it must be stressed that our finite channel is 
not designed to measure material parameters. There is 
a visible disagreement between the data and the fitted 
curves, which are based on an idealized theory. The solid 
lines are thus basically intended to guide the eye. 
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FIG. 9. Dispersion relation for the waves travelling 
counter-clockwise. 
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FIG. 10. Dispersion relation for the waves travelling clock- 
wise. 

In increasing the magnetic field the slope of the dis- 
persion relation is reduced. Because of the experimental 
difficulties described above, we have not attempted to 
obtain different wave numbers for one given driving fre- 
quency, but have rather plotted only the wave number kp 
of the largest peak. The prediction of a non-monotonic 
dispersion relation is thus not directly verified by the data 
points presented in the plot. The fit for B — 0.97 B c in 
Fig. |l0| shows a local minimum and makes the interpre- 
tation of the data points in terms of a non-monotonic 
dispersion relation very likely. Due to the inhomoge- 
neous magnetic field the corresponding measurement of 
the waves travelling counter-clockwise was not possible 
for this value of the magnetic field. 

FIG. 11. Space-time-plot for B = 0.96 B c , and T D = 2240 
Ths = 143.36 ms. 

In order to demonstrate that more than one wave num- 
ber is indeed present in the parameter range of interest, 
we finally present in Fig. |ll| a measurement for B = 0.96 
B c and a driving period of To = 143.36 ms. It clearly 
shows more than one wave number. In this case the rela- 
tively strong attenuation is caused by the longer driving 
period (compared to Fig. |J). The asymmetry between 
the left and the right hand side is believed to be caused by 
the inhomogeneity of the magnetic field, which becomes 
more dominant close to the critical value B c . The cor- 
responding spectrum is shown in Fig. |T^. The spectrum 
might indicate that we reach the limits of our measure- 
ment procedure here. It is fairly complex, presumably 
partly due to the effects of a nonlinear optical response, 
and effects of field inhomogeneities. In any case at least 
two peaks are clearly visible for the wave which travels 
clockwise. 
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FIG. 12. Spectrum obtained from Fig. [n]. Only wave num- 
bers corresponding to a frequency of T^ 1 — 6.98 Hz are 
shown. 
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IV. SUMMARY AND CONCLUSION 



By using a computer-controlled stroboscope algorithm 
we have measured non-monotonic dispersion relations for 
a magnetic fluid below the onset of a static pattern form- 
ing instability. Due to the finite geometry of our annular 
channel a quantitative comparison of the experimental 
results with a theoretical calculation is not possible. Our 
experiments clearly demonstrate a qualitative similarity 
with the theory obtained for a inviscous fluid in an infi- 
nite layer. Thus, the intriguing possibility arises to para- 
metrically excite PJlO[] three different wave numbers with 
one single driving frequency, and to study the nonlinear 
interaction of the ensuing surface waves. 
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